Layzer-Irvine equation: new perspectives and the role of interacting dark energy 
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We derive the Layzer-Irvine equation in the presence of a homogeneous (or quasi-homogeneous) 
dark energy component with an arbitrary equation of state. We extend the Layzer-Irvine equation 
to homogeneous and isotropic universes with an arbitrary number of dimensions and obtain the 
corresponding virial relation for sufficiently relaxed objects. We find analogous equations describing 
the dynamics of cosmic string loops and other p-branes of arbitrary dimensionality, discussing the 
corresponding relativistic and non-relativistic limits. Finally, we generalize the Layzer-Irvine equa- 
tion to account for a non-minimal interaction between dark matter and dark energy, discussing its 
practical use as a signature of such an interaction. 



I. INTRODUCTION 

The Newtonian energy conservation equation when 
generalized to an expanding cosmological background be- 
comes 



(1) 



where a dot represents a total derivative with respect to 
physical time, H is the Hubble parameter, E = K + U , 
with K and U being the peculiar kinetic and gravita- 
tional potential energies, respectively, of a system of non- 
relativistic particles interacting through gravity. This 
equation was derived independentely in [J 01 and it is 
known as the cosmic energy or Layzer-Irvine equation 
(see also Q)- Eq. ([!]) is valid throughout the entire pro- 
cess of structure formation and in the E — limit one 
recovers the usual virial relation K = —U/2 that holds 
for collapsed objects that have reached the state of hy- 
drostatic equilibirum. 

The Layzer-Irvine equation has been establishing itself 
as one of the most renowned equations of modern cosmol- 
ogy with its many applications including the determina- 
tion of the matter density, cluster mass and size, and 
the galaxy peculiar velocity field 043- More recently, 
some authors [8l-[Tl| have been using the Layzer-Irvine 
equation as a tool to detect a possible non-minimal in- 
teraction between the dark matter (DM) and the dark 
energy (DE) which, together, account for approximately 
96% of the energy content of the Universe today [l^, US] 
and whose fundamental nature is still largely unkown. 
The existence of such an interaction would in general 
invalidate the energy balance dictated by Eq. ([IJ. Con- 
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sequently, by measuring the properties of sufficiently re- 
laxed structures, such as galaxy clusters, one may expect 
to be able to detect a signature of an interaction between 
dark matter and dark energy through deviations from the 
usual virial relation 8- llJ. 

In this paper, our main goal is to present a broad dis- 
cussion of the Layzer-Irvine equation in a generalized 
framework (see also [11] )• In Sec. [ll] we start by de- 
riving the Layzer-Irvine equation in the presence of a 
homogeneous dark energy component, extending it to 
Friedmann-Robertson- Walker (FRW) cosmologies with 
more than three spatial dimensions. In Sec. Illll we show 
that the dynamics of cosmic string loops and other p- 
branes of arbitrary dimensionality are described by anal- 
ogous equations, discussing the corresponding relativistic 
and non-relativistic limits. In Sec. IIVI the Layzer-Irvine 
equation is generalized to the case where the dark mat- 
ter is non-minimally coupled to the dark energy back- 
ground and the implications of such coupling are dis- 
cussed. Finnally we conclude in Sec. |Vl 



II. NON-INTERACTING HOMOGENEOUS 
DARK ENERGY 



There is now overwhelming evidence for the Cosmolog- 
ical Principle which states that the Universe is homoge- 
neous and isotropic on cosmological scales. According to 
Birkhoff 's theorem, in the context of General Relativity, 
the gravitational field must vanish inside a spherical sym- 
metric shell, which is in agreement with the Newtonian 
result. This allows for the use of Newtonian mechanics 
in the study of the evolution of matter density fluctua- 
tions on scales much smaller than the Hubble radius. In 
this limit, the Lagragian for a system of point mass dark 
matter particles of mass m^, whose trajectories are given 
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by r.; = a(i)xi may be written as 



(2) 



where 



—rriivf 
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dt 
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(3) 

is the kinetic energy associated with the mass TOi, a is 
the scale factor, t is the physical time, Xj are comoving 
coordinates, Xi = |xi|, — Ti — Hvi is the peculiar 
velocity, H ~ a/ a is the Hubble parameter, Vi = jv^j and 
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^2^G(1 + 3^K^^^,^ (4) 



is the potential energy associated with the mass mi. 
Here, w = Pw/ Pw is the equation of state of the dark 
energy which we assume to be homogeneous (p^ and 
are respectively the dark energy density and pressure). 
We shall assume that the background evolution of the 
universe, given by a(t), is fixed, depending only on the 
average values of the density and pressure of dark matter 
and dark energy. Consequently, in Eq. ^ it is sufficient 
to consider only the inhomogeneous contribution to the 
gravitational potential energy due to DM-DM and DM- 
DE interactions. 

By performing the canonical transformation 



L -aa > \ 

dt \2 ^ 



the Lagrangian may be written as 
L = K-\J , 

with 

i 



where 



Ui = Ui + —adrriixf . 



(5) 

(6) 

(7) 
(8) 

(9) 



The first term on the right hand side of Eq. Q can 
be separated into two components: the contribution due 
to the homogeneous matter distribution plus the contri- 
bution due to the matter density perturbations. Then, 
using Eqs. dH) and ([9]), one obtains 



U = - — 



G 
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_Pm + P-ffi (1 + 3w)] + - 
[Pm(r) -pm] [Am(r') - A 



V (rriixf) 



^-dhd^r' 



(10) 



where pm is the matter energy density, pm is its aver- 
age value and the last equality was obtained using the 
Raychaudhuri equation given by 

^(p,„ + (l + 3u.)p^) = 0. (11) 



a 
a 



The Hamiltonian is given by 



Pi 

2m, 



with Pi 



rriiVi and the classical energy equation is 

dH _ dn 
H ^ 'dt' 



(12) 



(13) 



where the partial derivative with respect to time is com- 
puted at fixed particle comoving coordinates x^ and co- 
moving momenta Pi/a = niiXi. This way, one has 
U (X and K cx a~^. Consequently, using Eq. (fT3)) 
one finally obtains 



E + H{2K + U) = 0. 



(14) 



This shows that the minimally-coupled homogeneous 
dark energy does not explicitly enter the Layzer-Irvine 
equation. The effect of dark energy is felt only through 
the impact it has on the evolution of the Hubble parame- 
ter H. This generalizes the result in [l^ (where w = —1) 
to any homogenous dark energy form. For relaxed objects 
with E — one obtains the usual virial relation 



K = 



U 
2" 



(15) 



Extra dimensions 



It is interesting to generalize the above result to a A'' -|- 
1-dimensional FRW universe with N > 2. In that case, 
at fixed particle comoving coordinates x^ and comoving 
momenta Pi/a — miXi, one has U oc a~^"^^ and K cx 



which leads to 



E + H {2K + {N - 2)U) = 0. 



(16) 



Taking the case of sufficiently relaxed objects, for which 
_B = is a good approximation, then the virial relation 
becomes 



{N 



K 



which reduces to Eq. if TV = 3. 



(17) 



III. P-BRANE DYNAMICS 

The dynamics of maximally cosmic strings loops, do- 
main walls as well as higher dimensional p-branes in a 
cosmological background has been studied in detail in 
\vf\ . This work has recently been extended to ac- 
count for the dynamics of cosmological p-brane networks 

MM- 
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A. Cosmic Strings 

In the absence of non-gravitational interactions (as well 
as gravitational radiation backreaction) the evolution of 
the total energy i? of a cosmic string loop is given by [T6j 



with 



E = 2HE { ^ - v^ 



with 



E = fj,a J 'yds . 

-2 ^ Jv'^lds 
lids 



(18) 

(19) 
(20) 



where ^ is the energy per unit length, ds is the infinites- 
imal comoving arclength, v is the loop velocity at a par- 
ticular point and 7 = (1 — u^)^^/^. For very small loops 
(with E/^ ^ H~^) it is in general a good approximation 
to consider that the expansion has, on average, no impact 
on the total energy. Hence, the average over a sufficiently 
long time of the total energy {E)t and root-mean-square 
(RMS) velocity (5^)f = 1/2 is approximately constant. 

In the case of a very large non-relativistic loop, the 
total energy can be decomposed into the potential en- 
ergy associated with the loop length U = ijlL oc a, 
and the kinetic energy associated to the loop motion 
K = fiLv^ /2 (X {v oc a~^), where L is the physi- 
cal length of the loop. As a result, using Eq. one 
obtains. 



E + H{3K -U) = 0. 



(21) 



Eq. (|2ip is very similar to the Layzer-Irvine equation: 
in both equations the derivative with respect to physical 
time of the total energy is proportional to the Hubble 
parameter times specific linear combinations of the ki- 
netic and potential energy terms. Note however that in 
the case of non-relativistic cosmic strings one cannot set 
E = and therefore there is no analogy with the grav- 
itational virial relation. This happens because in the 
non-relativistic regime K ^ U so that E ^ U oc a. An- 
other diference is that, contrary to Eq. (j21l) that has 
a relativistic version (Eq. (IT5|) '). there is no relativistic 
generalization of the Layzer-Irvine equation. 



B. p-branes 

Analogously to the case of the Layzer-Irvine equation, 
we can also generalize the cosmic string case to higher 
dimensions. In the case of a p-brane, Eq. ([T5| generalizes 

to 



E = {p+ l)HE 



p+1 



(22) 



E = apQP J 'fdA, 



J v^-/dA 
JjdA ' 



(23) 
(24) 



where ap is the energy per unit p-dimensional area, dA is 
the infinitesimal comoving p-dimensional area, v is the p- 
brane velocity at a particular point and 7 = (1 — v^)~^^^. 
For very small p-branes (with {E/apY/^ <C -ff^^) the 
expansion has in general a very small impact on the 
time average of the total energy {E)t and RMS veloc- 
ity (w^)t — p/{p + 1) which are therefore roughly con- 
stant. On the other hand, for very large non-relativistic 
p-branes one has U = (TpA oc and K = cFpAv^ /2 oc 
U-2-P (^y a^p-i)^ where A — a^A is the physical p- 
dimensional area of the p-branes. As a result, the energy 
equation becomes 



E + H i{2+p)K -pU) ^0. 



(25) 



This equation generalizes Eq. (|2ip to p-branes of arbi- 
trary dimension in N-|-l-dimcnsional homogeneous and 
isotropic FRW universes (with p < N). The similarities 
with the Layzer-Irvine equation are again very evident. 



IV. INTERACTING DARK ENERGY 

One of the ways to better understand the physics of 
dark energy is through its influence on the formation of 
large-scale structures in the Universe. In an accelerated 
Universe the characteristic timescale for linear perturba- 
tion growth may become large compared to the Hubble 
time. However, if dark energy and dark matter interact 
non-minimally, then dark energy influences the process 
of structure formation in a more active way, not only 
through its impact on the acceleration of the Universe. 

The coupling between dark matter and dark energy 
adds new source terms to the usual Layzer-Irvine equa- 
tion 

[Mil (Eq. ©). These extra terms can be written, 
with all generality, as 



dK 

'at 



dU 



mHU, (26) 



so that the generalized Layzer-Irvine equation becomes 



E + H{{2- a)K -f (1 - I3)U) = . 



(27) 



The functions a{t) and I3{t) depend on the details of the 
process of energy and momentum transfer between dark 
matter and dark energy and are therefore model depen- 
dent [13]. For example, if dark energy decays into dark 
matter, then new particles with non-vanishing momen- 
tum may be continuously added to the system. This way 
it would be crucial for the computation of a{t) and P[t) 
to know not only the rate of energy transfer but also 
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the initial RMS velocities of the new particles. On the 
other hand, the coupling might also occur through the 
dependence of the mass of the dark matter particles on 
the value of the dark energy field (see, for example, [2l|). 
In this paper, the model dependence associated with dif- 
ferent choices of coupling models is incorporated in the 
freedom to choose the evolution of the parameters a{t) 
and 13 (t). 

In Q the case with a = was considered, with 
/? being related to the coupling strenght. It was argued 
that P could be determined by measuring the kinetic and 
potential energy of sufficiently relaxed structures such as 
galaxy clusters. The homogeneous dark energy case with 
a — j3/2 has also been considered in If a and /3 

are constants then the virial relation obtained assuming 
hydrostatic equilibrium [E = 0) is given by 



2-a ' ' ' 

We note however that in the presence of such an interac- 
tion one cannot, in general, assume hydrostatic equilib- 
rium. This can only happen if a and /3 are constant or, 
according to Eq. (P7| , if their evolution is given by 

m = {E+{l~a{t))K)/U , (29) 

with constant E, which does not happen in general. 
These two cases are very special and consequently, in the 
presence of an interaction between dark matter and dark 
energy, gravitationally bound systems are not expected 
to reach virial equilibrium. As a result, deviation from 
the usual virial relation in galaxy clusters is therefore 
a general signature of a non-minimal coupling between 
dark matter and dark energy. 



V. CONCLUSIONS 

In this paper we studied the Layzer-Irvine equation 
and discussed some of its generalizations. In particular, 
we derived the Layzer-Irvine equation in the presence of 
a general homogeneous dark energy background showing 
that the final form of the equation is not affected explic- 
itly by the dark energy component. We further gener- 
alized the equation and the virial relation to FRW cos- 
mologies with -I- 1 dimensions (with N > 2). We have 
also demonstrated that the macroscopic dynamical en- 
ergy equations of cosmic string loops and other p-branes 
of arbitrary dimensionality are, in the non-relativistic 
limit, analogous to the Layzer-Irvine equation. Finally, 
we generalized the Layzer-Irvine equation to account for 
a non-minimal interaction between dark matter and a 
homogeneous dark energy form. We have shown that, in 
general, gravitationally bound systems are not expected 
to reach hydrostatic equilibrium in the presence of a cou- 
pling between these two components. This constrasts 
with the usual assumption in the literature where the 
equilibrium relation E = is assumed a priori. Hence, 
a non-minimal coupling between dark matter and dark 
energy will generally lead to the breakdown of the usual 
virial relation K = —U/2, providing a crucial signature 
of such an interaction. 
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